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ABSTRACT 


Computer — aided control system design, using quantitative 
feedback theory (QFT), for continuous time, linear time invariant, 
single-input single-output, systems with parametric uncertainty, 
has been tackled in this thesis. The general parametric uncertain 
model assumed here involves a transfer function with its 
coefficients given as polynomials over a set of system parameters, 
varying in a real hyper-rectangle. A modified QFT algorithm, 
based on the extremal, rather than, as is done traditionally, the 
entire, system frequency response, has been implemented here. The 
attendant frequency response extrema computation is based, for the 
general polynomial uncertainty case, on an uniform grid search 
scheme. For a special case of systems with polynomial 
uncertainty, viz. the linear affine uncertainty systems, the 
extremisation is based on a frequency response mapping theorem. 
This results in a computationally efficient algorithm, relative to 
the grid method, for the linear affine uncertainty case. The 
consequent improvement in the computational aspects and the ease 
of applicability of the basic QFT design paradigm, achieved using 
the modified QFT algorithm, relative to the traditional QFT 
approach, is demonstrated through the solving of five design 


examples. 
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CHAPTER 1 
INTRODUCTION 


1«1 The Feedback Issue 

Technological systems like electronic amplifiers, process 
plants, air and space craft, ships, etc., are designed to perform 
in a certain required manner. But it is usually observed that 
these systems exhibit a behaviour that fluctuates about a desired 
or nominal behaviour. The reasons for this can essentially be 
traced to the following factors - 

(a) Mathematical model uncertainty - The design process 
involves the use of mathematical models which for 
reasons of ease of understanding and applicability, are 
just approximations to the actual system. .The resulting 
loss of system information renders the system behaviour 
about the designed for value unpredictable or uncertain. 

(b) External factors - The design process also assumes a 
certain specific set of inputs to the system. However, 
external factors (also known as noise or disturbances) 
like e.g. electromagnetic interference, ambient 
temperature fluctuations, wind speeds, etc.,? cause the 
inputs to change randomly about the set values thereby 
inducing uncertainty in the system output. 

Certainly remedial measures would be called for if the system 
behaviour uncertainty (in the form of output fluctuations) is "too 
much" i.e. cannot be tolerated. Feedback control is a technique 
which specifically and effectively addresses this problem viz. 
that of reduction of system behaviour uncertainty in the face of 
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model and input uncertainty CH0R63]. Feedback essentially 
involves monitoring the time flow of the system characteristic to 
be controlled like e.g. output of a voltage supply, temperature an 
a chemical reactor, motion of an aircraft or ship, etc., and based 
on its behaviour rectify the input responsible for this system 
characteristic like e.g. input voltage, fuel flow, etc. 
Monitoring and rectification could involve processes like 
amplification, transduction, etc. . Setting up these monitoring 
and rectification sub-systems is what control system design is all 
about . 

1.2 The Feedback Design Problem 

A commonly used block diagram representation of the feedback 
idea is as shown below - 



Figure 1.1 The basic feedback configuration 

Here, P is the given system (also known as the plant) and "r" 
IS the desired output. The actual system output "y" is compared 
with "r" and the resulting error "e" ( = r - y) is used to drive 
a corrective element G (also known as a compensator) which in turn 
drives P in such a manner as to reduce the error "e" to zero. 
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"dl" and "d2" are external plant input and output disturbances, 
respectively . 

Synthesis of a "low cost" G constitutes the feedback design 
problem. The cost of G is specified in terms of the complexity of 
G with regard to its gain and bandwidth. Given a specific problem 
i.e. a plant P and a set of specifications characterising the 
output "y", many approaches are available for synthesising G. 
However, we could possibly classify all the existing design 
approaches into the following two broad classes — 

CLASS 1 Classical approach — This is the approach 

encountered in most control engineering texts like 
e.g. CBOWbS], CTAK663 , |:H0R63], CH0SB9], [D AZ60] , 

to name a few. Here, design is done for only 

a single model of P (known as the nominal model). 
In this method, pioneered essentially by Hendrik 
Bode CB0D45], uncertainty is indirectly and crudely 
accounted for through the specification of two 
design parameters, viz. the gain and phase margins 
CB0D45], CD'AZ60]. The specification of these 
margins is, however, highly heuristic without any 
rigorous analyti"al basis. Thus improper 

specifications are possible which could lead to a 
great deal of trial and error (commonly referred to 
as "tuning the controller") in the design process. 
This aspect makes the following approach more 
attractive. 

CLASS 2 Modern approach - Here, the model uncertainty, 
quantified in an appropriate manner, is explicitly 
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utilised. This feature essentially eliminates a 
great deal of cut and try that may have to be 
resorted to, when alternatively, the class 1 
approach is used. Mathematically rigourous and 

quite reliable theories of LQG— LTR CD0Y81], H 

oo 

control CFRA87] , analysis / synthesis CSTE91] and 
quantitative feedback theory (QFT) CH0R91] fall 
under this class. 

QFT, a frequency domain technique, pioneered by Isaac 
Horowitz, differs from the other techniques in the way uncertainty 
is quantified. In the frequency domain, there are mainly two 
types of model uncertainty, viz. structured and unstructured. The 
uncertainty in the system induced by variable system parameters 
like e.g. variable value passive electrical components, variable 
value gain factors, variable value inertia elements, etc., is said 
to be structured . This uncertainty makes itself felt essentially 
in the low frequency range. The uncertainty induced by parasitic 
elements and unmodelled dynamics, which make themselves felt at 
high frequencies, is said to be unstructured . References in this 
regard are CD0Y81], CKWA93], CH0R91]. QFT utilises the structured 
uncertainty representation as aga nst the use of the unstructured 
representation by the other techniques. 

As a design technique QFT is quite appealing relative to the 
other techniques on account of the following reasons - 

(1) The structured uncertainty representation is more 
realistic relative to the unstructured representation 
CBHA87]. Also the design effort (in terms of the gain 
and bandwidth of G) is usually exercised over the low 
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frequency range, over which structured uncertainty is 
more predominant than the unstructured representation. 
These factors and the explicit usage of structured 
uncertainty make QFT well suited to tackle realistic 
design problems. Though it is possible to handle low 
frequency structured uncertainty using the other 
methods, this however, involves conversion of the 
structured representation into unstructured CLAU86] and 
in the process could result in overdesign i.e. 
unnecessary design, or worse still, fail to give a 
solution for an otherwise QFT solvable problem CBAI91]. 

(2) In the multi-input, multi-output (MlhO) QFT problem, the 
issue of guaranteeing specifications for the individual 
transfer function entries of the closed loop system 
transfer matrix is directly addressed to. This feature 
IS absent in the other techniques CH0R82]. 

(3) In the other techniques, synthesis of G is very much 
dependent on the choice of a nominal model for the 
system. This necessiates a search for an optimal 
nominal model which wv.ald result in the cheapest 
controller G. In contrast, QFT is quite insensitive to 
the choice of a nominal model for the system P. In fact 
as shown in this thesis, a nominal model need not be 
chosen at all for the QFT design procedure. 

(4) The trade-offs involved between the performance 
specifications and the cost of G are more transparent 
in QFT than in the other techniques. 



1.3 S^ope of t.hie Thesis 

Software for handling QFT design for single— input , 
single-output (SISO) continuous time, linear time invariant (LTI) 
systems has been developed in this thesis. 

The attempt here was to fulfil the need felt for a computer 
aided QFT design scheme. To date, no commercial software for the 
same is available CSPEC93]. The reliability of this software was 
tested by using it to solve five examples. These five examples 
were so chosen as to cover, as far as possible, a broad spectrum 
of SISO systems and design problems encountered in real life. 
They include an unstable system and a nonminimum phase system as 
well. This software, when used in conjunction with MATLAB 

(including the CONTROL TOOL BOX) CMATS7], should provide a very 
easy to use tool for tackling real life feedback system design 
problems as well as assist in teaching QFT. Further research in 
QFT can also be facilitated using this software. 

1.4 Organisation of the thesis report 

This thesis report is organised in the following manner - 

Chapter 2 gives a formal mathematical definition of the 
feedback design problem tackled by QFT. 

Classical QFT as pioneered by Marcel Sidi and Isaac Horowitz 
CHQR72], IS described in chapter 3. 

Chapter 4 contains a description of the algorithm used in 
implementing the computer — aided QFT design scheme of this thesis. 

The solution for the five examples obtained using this 
software are presented in chapter 5. 

Finally, a summary of the work done and scope for future work 



are discussed in chapter 6. 


Appendix A contains a description of the software and its use 


explained by means of an example. 
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CHAPTER 2 

QFT PROBLEM STATEMENT 

In this chapter we will give a formal mathematical 
description of the feedback control system design problem handled 
in this thesis using QFT. 

2. 1 System Description 

QFT utilises a model for the plant P incorporating the 
structured uncertainty representation. As mentioned in section 
1.2, chapter 1, this type of uncertainty arises essentially due to 
the fluctuations in system parameter values about their mean or 
nominal values. More specifically, these system parameter 
variations give rise to a family of transfer functions given as - 
P(s) = -C P(s,q) , q € Q) 2.1(a) 

where 

Q={q€lR’"|q"<q <q^,l<i<n} 2.1(b) 

where q is an n-tuple of real uncertain system parameters. Each 
transfer function P(s,q) belonging to this family can be 
cons dered to have a form given as - 

Tn 

Z bL(q) s"^’- 

P(s,q) = (2.2) 

n 

E ai(q) s''"’’ 

1=0 

Thus all the P(s,q) belonging to P(s) can be considered to be 
correlated to each other via the functional dependence of the 
coefficients at and bt on the uncertain parameter vector q. 


A set of uncorrelated transfer functions obtained by 
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linearisation of a nonlinear system model about several operating 
points, IS also sometimes considered to constitute a structured 
uncertainty model for the system involved [AZV92]. However, in 
this thesis we were not concerned with this type of structured 
uncertainty models. Consequently, it should be noted that the 
software developed in this thesis cannot handle this type of 
design problems. 

2.2 Design Specifications 

Feedback control systems can essentially be classified into 
the following two types CTAK66] - 

(1) Tracking control systems - Here the contolled output 
y(t) (refer figure 1.1) is required to follow or track a randomly 
varying reference input r(t). Here, the most desirable system 
property, obviously, is that y(t) should follow the variations in 
r(t) as rapidly and accurately as possible. Generally used 
performance measures for this property are - 

(i) Rise and settling times of the system unit step 
response . 

(ii) Peak overshoot of the system unit step response. 

A generally used single loop tracking system is as shown 
below - 



Figure 2.1 The general feedback configuration 
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"I'his system is the same as that given in figure 1-lj except 
for the additional feature of a pre— filter F(s) . The compensator 
G(s) reduces the variation in the output y(t) due to system model 
and inputs uncertainty, while the pre-filter F(s) helps in 
shaping y(t) to satisfy the tracking specifications. 


(2) 

Regulators - 

Here the output 

y( t) 

is sought 

to 

be 

maintained at a fixed 

value throughout 

the 

operation 

of 

the 

system . 

Consequently , 

the reference r(t) 

IS a 

constant . 

Here, 


when the system is subjected to a disturbance, the controlled 
ouput y(t) should retrieve its original value as quickly as 
possible and with the least steady state error possible. Also, 
the resulting transients should be sufficiently damped. Generally 
used performance measures are - 

(i) Settling time and peak overshoot of the system 
impulse response. 

(ii) Steady state error of the system unit step 

response . 

Uncertainty in the system model and inputs, forces us to 
define a tolerance band for the values assumed by the performance 
measures mentioned above. For example, for a unit step response 
we may desire a specification such that 
Peak overshoot Mp e [Mp , Mp] , 

Rise time tr e [tr , tr], etc. . 

Since QFT IS a frequency domain technique, these performance 
specif ications , which are in the time domain, have to be first 
translated to the frequency domain. Unfortunately, however, to 
date, there is no rigourous technique to carry out this 


translation. Consequently, 


a more or less heuristic technique 
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IS adopted -for this purpose CAZV92], CD AZ88], CASH823. 

Here, the key lies in the excellent correlation between the 
time domain and frequency domain responses of a second order 
system and the fact that a higher order system s response can be 
approximated by that of a second order system s CTAK66] , 
[D AZ883. 

Two transfer functions, Tu(s) and Tl(s) , are found which 
approximately define the upper and lower bounds, respectively, of 
the given time domain performance measures. If the closed loop 
transfer function from r(t) to y(t), viz. Tyr(s), is such that 

0 < I Tl(jw) I < I Tyr(j<o) I < I Tu(j<o) |, V « > 0 (2.3) 

then, as a rule of thumb, the time domain performance measures 
will be satisfied. It should however be noted that this may not 
be always true. Consequently, some amount of cut and try could be 
required in effecting the required time domain to frequency domain 
conversion. To facilitate this trial and error process, two 
second order transfer functions are chosen as an initial guess for 
Tu(s) and TKs) CHOR723, CAZV923, CD-AZ883, CASH82D. Simple 
poles and zeros can then be added to the initial guess to 
represent the time domain specifications more closely. This 
process can be understood better via the following example - 

EXAMPLE 2.1 - Suppose that we have the following time domain 
specifications (refer example 5.1, chapter 5) : 

SPl ! Peak overshoot Mp of the closed loop unit step 
response (refer figure 2.1) should be less than 
20 7 .. 

SP2 : 27. settling time ts of the closed loop unit step 
response should be between 6 and 8 seconds. 
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For the time domain to frequency domain translation, 
an initial guess for Tyr(s) and Tyr(s), the "boundary" functions 
here, was taken to be a second order transfer function of the form 
given below — 

Tguess(s) = 


2 

<un 


2 2 
S + (2^0>n)s + con 


where ? = damping factor and <A>n = undamped natural frequency. 
Using the formulae, 

- 


Peak overshoot (V.) = e 




Settling time % 


4.0 


if con 


and MATLAB , the following two boundary models were obtained. 


T yr ( s ) = 


1.08 


s + 0.94s + 1.08 


,with Mp = 20% and ts = 8 sec 


Tyr(S) = 


1 .8 


s + 1.248s +1.8 


,with Mp = 18% and ts = 6.2 sec 


such that 


1 Tyr(jco) ( < | Tyr(jco) | , Vco > 0. 


To increase the spread etween | Tyr(jco) j and | Tyr(jco) | with 
increasing frequency, without really compromising the time 
responses, for reasons made clear in the next chapter (see section 
3.3, chapter 3), we added a pole to Tyr(s) and a zero to Tyr(s), 
resulting in 

, , (0.216s + 1.08) 

Tyr(j<0) = 


s + 0.94s + 1.08 


with Mp = 21% and ts = 7.8 sec 


and 
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was obtained through trial and error aided by MATLAB such that 
I TyrCjco) I < | Ttest(ja>)| < | Tyr(jco)|, Vw>0. 

The test transfer function, as expected, had a Mp = 1127 . and a 
settling time of 6.75 seconds. Consequently, we can now assume 
that if the closed loop transfer function Tyr(s) is such that it 
satisfies equation 2.5(a), with the above given Tyr(s) and Tyr(s) 
functions then the time domain specifications of SPl and SP2 will 
be met. This assumption is more or less vindicated in examples 
5.1 and 5.2, chapter 5. 

Similar frequency domain specifications of the form given in 
equation 2.3, but with only the upper bound Tu(s) defined, can be 
used to specify the system disturbance response [D AZ88], CAZV923. 


2.2.1 A Digression 

In the QFT literature, e.g. CH0R72D, CAZV92], one can 
find attempts to satisfy time domain specifications of the 
form - 

l(t) < y(t) < u(t), Vt>0, VqeQ (2.4) 
where l(t) and u(t) define time domain bounds on the 
desirable system output y(t). But as pointed out in CJAY933 , 
CBAR93] , this is quite a tough problem to solve as conversion 
of the equation 2.4 - type specifications to that in equation 
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2 3 IS quite difficult and highly heuristic with not much 
theoretical support The methods most commonly suggested are 
of the sort described in example 2 1 above CASHS2], CH0R72] 
Also, CASH82] reports a Monte Carlo determination of the 
allowable frequency response spread corresponding to equation 
2 4, given the knowledge of the model parameters and their 
variances. CKRI77] is another related reference However, 

the use of these methods is not widespread 

But note that equation 24 is essentially a tracking 
specification Consequently, we feel that it would be much 
more easier to first convert the equation 24 - type 

specifications into an equivalent set of specifications on 
the unit step response of the system, e g in terms of bounds 
on the rise time, peak overshoot, settling time, etc 
Subsequently, these step response specifications can then 
be converted into the equation 23 - type specifications, 

relatively easily, as illustrated in example 2 1 above 
However, this type of a design problem has not been 
considered in this thesis 

2 3 Problem Statement 

Summarising the above given discussion, we can now formally 
state the QFT design problem tackled in this thesis, as follows - 
Consider an uncertain continuous time system, with a linear 
time invariant model of the form given in equation 2 2 and 
embedded in a feedback configuration of the form given in figure 
2 1 The QFT design problem, then, is to obtain fixed (i e with 
no uncertainty), physically realisable, LTI compensator G(s) and 
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pre— filter F(s) such that the resulting closed loop system shown 
in figure 21 is stable and the following frequency domain 
specifications are satisfied - 

I TyrCjto) I < I Tyr(joj) j < | Tyr(jt»j) | 2 5(a) 

I Tydi(jco) I < I TydiCjto) I 2 5(b) 

I Ty<i 2 (joj) I < I Ty<i 2 (jw) 1 2 5(c) 

V oj > 0 , 


where 


Tyr ( s ) 


y(s) 
r (s ) 


F(s)L(s) 
10+ L(s) 


2 6(a) 


T y<bl ( S ) 


y(s) 

dl(s) 


P(s) 

10+ L(s) 


2 6(b) 


Tyci 2 (s) d2(s) ~ 10 + L(s) ^ ^ 

where 

L(s) = P(s)G(s) IS the (open) loop gain 


A description of the application of the OFT design technique 
IS what follows in the next chapter 


16 


CHAPTER 3 

THE QFT TECHNIQUE 

In this chapter we will describe the QFT design procedure 
pioneered by Marcel Sidi and Isaac Horowitz [H0R72] 

3 1 Essence of the QFT design procedure 

Uncertainty in the system model results in a fuzzy phase 
magnitude Bode plot More specifically, at each frequency wo the 
frequency response of the system consists of a set of points in 
the complex plane, given as - 

TmP(too) = { P(jwoq), q e Q> (3 1) 

as shown below - 



Figure 3 1 Uncertainty template at wo 


This set IS known as a template CH0R723 or a value set CFU90] 
of the system at the frequency ojo 

Now, the closed loop transfer function from the reference 


input r(t) to the output y(t) is (refer figure 2.1) given as - 

F(s)L(s) 


Tyr(s) = 


10+ L(s) 


(3 2) 


where 


L(s) 


P(S)G(S) 


IS the (open) loop gain 



17 


Using the log magnitude notation, where 

LmCT(ja))] = 2(Z)log(| T(jco)|) (3 3) 

with "log" denoting the common logarithm i e to the 
base 10, we get 

LmCTyr(ja>)3 = LmCF(j<o)3 + LmCL(j<o)/(l 0 + L(j<o))] 
with LmCL(j<*>)] = LmCP(j<«j)] + LmCG(jaj)] (3 4) 

Uncertainty induced variation in causes variation in 

Tyr(j<o) such that - 

ALmCTyr( j to) ] = ALmCL ( j to) / ( 1 0 + L(jto))] 
with iiLmCL(jto)] = ALmCP(joJ)] (3 5) 

Note that G(s) and F(s) are assumed to be fixed (refer 
section 2 3, chapter 2) 

From equations 3 1-3 5 it is quite clear that in the Nichols 
chart CD AZGS] the loop gain template 

TmL(too) = CP( jtoo,q)G( jtoo) , q e Q> (3 6 ) 

has the same shape as the plant template TmP(coo), but is 
translated vertically by LmCG(J<*>o)] and horizontally by 
argC6( j coo) ] , where argCB(jcoo)] is the phase of G(j<oo), as shown 
below - 



Figure 3 2 Template movement in the Nichols chart 
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Also, observe that it is in terms of the magnitude deviation 
of the complementary sensitivity function 

Tc*(s)=L(s)/(10+L(s)) (37) 

that the closed loop magnitude deviation is expressed (see 
equation 3 5) Consequently, given a loop gain template TmL(coo) 
in the Nichols chart, the corresponding magnitude variation in the 
closed loop gain at the frequency a)o can be directly read off the 
chart via the constant M - contours This aspect and the feature 
of permanency of the loop gain template shape irrespective of the 
value of G(j<»>o) are made use of in the QFT algorithm to be 
discussed in the course of the following discussion 
Now, the specifications 2 5(a) dictate that 

ALm[Tyr( j «) ]meut i Lm[ Tyr ( j <i>) ] - LmC Tyr ( j o>) ] , V(0 S 0 (3 S) 
From equation 32 it is quite clear that if 
(10 + L(jcoo)) ?^0, Vq e Q, then 

Tyr(jcoo) > F(jtOo) as | G(jtoo) | t Oo 

ALmCTyr( J Wo) ] ► 0 

In other words, ALmCTyr( j <oo) ] can be made arbitrarily small 
by choosing G(jcoo) with LmCG(jwo)3 sufficiently large This 

essentially means that there exists a Gi|((jo>o) value such that 
equation 3 8 will get satisfied for all G(jwo) values having the 
same phase as G*(jwo) but satisfying the inequality 
LmCG(jo)o)] > LmCG*(jwo)] 

A formal proof of this statement is given in CBAI88] 

Such G(jwo) values will be said to be admissible performance 
Mise The design objective is to synthesise such a performance 
wise admissible compensator G(s) which additionally ensures closed 
loop stability (i e. is admissible stability wise). 
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This G(s), however, only guarantees satisfaction of equation 
3 8 and does not ensure satisfaction of the tracking 
specifications of equation 2 5(a) This gets done by the 

pre-filter F(s) For this the pre-filter magnitude at uo should 
be such that 

LmCTyr( J <ao) ] - LmC Tyr( J coo) ] < LmCF(ja>o)] 

mtn 

£ Lm[Tyr( J<0o) ] - LmCTyr( J <oo) 3 (3.9) 

mou< 

Any minimum phase transfer function F(s) which satisfies 
equation 3 9 for all uo > 0, would provide a solution for the 
pre-filter F(s) 

On similar lines, a performance wise and stability wise G(s) 
IS also sought for the disturbance specifications of equations 
2 5(b) and 2 5(c) Note that a pre-filter is not needed for 
ensuring these specifications 

The design procedure, at the very outset, thus requires a set 
of frequencies, O, to be identified such that 

(10 + L(ja>))?<0, Va> e O 

This set IS known as the trial frequencies set For 

computational solvability we naturally consider a finite 
Though the selection of fl can be arbitrary, a common practice is 
to include frequencies lying in the DC to <oh range, where «h is a 
high frequency value such that the structured uncertainty model 
P(s) (refer equation 2 2) tends asymptotically to k/s" 

(k variable) for all frequencies <o > ojh [D AZ883 , CAZV92] The 

mam reason for this choice is that control effort is usually 


concentrated in this low frequency range 


zo 


3. 2 The Sxdj. — Horowit.z QFT algorilhin 

This algorithm has the following steps - 

STEP 1 Setting up the frequency domain performance 

specifications - This essentially involves a 
a conversion of time domain specifications to the 
frequency domain, essentially on the basis of a 
semi-heuristic procedure of the kind described in 
example 2 1, section 2 2 chapter 2 
STEP 2 Specification of bounds for a performance admissible 
G(s) - At the very outset, an appropriate finite 
set of trial frequencies O is selected Then, at each 
trial frequency uo at a fixed assumed phase of G(j{*to), 
VIZ 0(coo), a value G*(jwo), with phase <^(coo) , is found 
such that for all G(jwo) values with phase 0((oo) and 
I G(j6jo) I > I G*(jcoo) |, equation 3 8 or 2 5(b) or 2 5(c) 
(according to the problem type), is satisfied Several 
such "bounding" G*(j<oo) values are found by varying 
^(coo) over some desired phase interval of length 360 
degrees (e g the interval C-360'*,0®3) 

As it IS the loop gain which is important in the 
design procedure, these G*(jooo) values are consequently 
used to define equivalent lower bounds 

L*(ja)o) = P( jtoo q)G» ( jojo) , on a nominal loop gain value 
Lo(jo>o) defined, arbitrarily, for some nominal parameter 
vector value q© e Q Thus, we get several such 

L#(j<oo) points which can be joined to form a closed 
continuous curve in the complex plane known as a B(<j 3) 
bound for the frequency <oo as shown below - 
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It follows quite easily that for G(s) to be 
admissible performance wise, the nominal loop gain 
Lo(j wo) should lie above i e outside the corresponding 
B(wo) bound at each trial frequency wo, as shown in the 
Nichols chart below - 



If simultaneous specification of two or all 
equations 2 5 (section 2 3, chapter 2) is done then for 
every trial frequency o>o, we will correspondingly get 
two or three B(wo) bounds Then we, obviously, take 
the worst combination of these bounds to get a 
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composite B(<oo) bound, as shown below - 



The Nichols chart is very convenient in finding 
these B(a>) bounds For instance, at each assumed value 
of the phase of G(j<a>o), viz <^(a>o), the L#(juo) value is 
found by simply moving the corresponding loop gain 

template vertically up or down as need be, till, that 

position IS attained where the performance 
specifications just about get satisfied CH0R72] , 
CD AZ88] The amount of translation (in terms of dB ) 
done in attaining this position will then define the 
GiK(j(oo) value 

STEP 3 Loop shaping for stability - In this step, the 
stability of the closed loop system with the G(s) 
obtained from step 2 is tested by subjecting the 
resulting fuzzy open loop polar plot to the Nyquist 
encirclement test If the test is passed then go to 

step 4 , else modify G(s) to ensure closed loop 

stability and go back to step 2 to test for performance 
admissibility of the modified G(s) 
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This step involves the usual cut and try of the 
classical loop shaping method CAZV92] is a very useful 
reference in this regard 

STEP 4 Pre-filter design - If the pre-filter F(s) is not 
required then end the design procedure, else synthesise 
a minimum phase transfer function F(s) satisfying 
equation 3 9 This F(s) will then constitute a solution 
for the pre-filter and the design procedure can then be 
ended 

Thus QFT IS essentially an extension of the traditional loop 
shaping methodology to design problems with large parametric 
uncertainty 

It will be useful to note that it has been shown [GER80] 
that an optimal solution for G(s) is the one in which all the 
Ltf>(jcoo) points lie on their respective B(a>o) bounds This 

solution IS optimal in the sense that the resulting | L(jto) | will 
have the fastest possible roll-off rate (as the frequency 
increases) and in the process, the least sensor noise problem 
CH0R63], CH0R91] 

3. 3 Further issues in the design procedure 

3 3 1 High frequency behaviour of the B(a>) bounds 

The generation of the B(o>) bounds, for each trial 
frequency too, is done by varying ^(too) over an interval of 
length 360® This naturally creates the impression that the 
B(to) bounds for all the trial frequencies, encircle the 
origin, as shown in figure 3.3 That this may not be 
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necessarily true (and that an fact it is not desirable) can 
be explained as follows 

For most physical systems, PCs) ► 0, as 

|s| »■ 00 Also, as 6(s) is assumed to be proper or 

strictly proper (the assumption of real isabi 1 i ty ) we have 

L(s) ► 0, as |s| > 00 Consequently, we have 

Tyr(s) ► F(s)L(s) , as |s| *■ co 

^ ^mCTyr(joo)] ► 4LmCP( j<i>) ] , as o> »• oo 

This means that as a necessary condition for the 
satisfaction of the tracking specifications in equation 
2 5(a), we must have 

ALmCP ( J <*>) ] — Lm [Tyr ( j w) ] — LmCTyr( joo) ] , Vco i to# (3 10) 

where oo# is a very high frequency 

In other words feedback is really not needed for 
frequencies beyond w# Now, equation 3 10 can be shown to 

guarantee that as co ► oo , the BCco) bounds start shrinking 

and beyond co# cease to encircle the origin but persist in 
encircling the -1+J0 point in the complex plane as shown 
below CH0R733, CH0R76], CH0R79] - 



Figure 3 6 High frequency behaviour of the B(co) bounds 
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Consequently, we can define a closed contour L) enclosing 
the -1+J0 point but not encircling the origin, as shown in 
the above given figure, such that for all frequencies <o > to* 
the corresponding B(aj) bounds lie entirely inside the closed 
contour U This contour is known as the universal high 
frequency (UHF) boundary CH0R72] , [D AZ88], as avoidance of 
U for all oi 03$ ensures avoidance of the B(w) bounds for all 
the frequencies above 

The UHF boundary U has an additional significance in 
that for certain disturbance attenuation problems, a 
constraint on the damping factor is enforced via an M-contour 
constraint on the complementary sensitivity function Tcs(s) 
[H0R72] , [D AZaS] That is, the performance requirement is 
to satisfy 

I Tcs(jo>) I < M , Vw > 0, Vq e Q. 

This constraint can be enforced by defining the 
concerned M-circle as the UHF boundary 

The undesirability of persistent encirclement of the 
origin by the B(a)) bounds (i e for all co) can be appreciated 
due to the fact that it could then imply an unrealisable 
G(s) To see this, note that persistent B(o)) encirclement of 
the origin implies that the nominal loop gain Lo(s) cannot go 
to zero, asymptotically or otherwise If P(s,qo) is strictly 
proper then, to ensure that Lo(s) does not go to zero as 
|s| ► 00 , G(s) will have to be improper and hence 


unrealisable 
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332 Variable relative degree of the structured uncertainty 
model 


The 

relative degree ' e" 

of 

P(s) (refer 

equation 

2 2) IS the 

difference between 

the 

denominator 

and the 


numerator degrees i e e = n - m 

A variable relative degree implies the existence of some 
"disappearing" poles and zeros giving rise to a 
representation of P(s,q) of the form - 

n ( bis + 1 ) 

P(s,q) = Pi(s,q) (3 11) 

pj (ats + 1) 

L 

where Pi(s) has a fixed relative degree and each bt and a^ 
can take the value 0, thereby giving rise to the phenomenon 
of disappearing poles and zeros and the consequent variation 
in the relative degree e CH0R793 

Now, as P(s,q) *■ k/s* (k variable), for|s| * oo , 

if e(max) > 4 , then the loop gain template could have a 
phase variation of more than 360*^ and consequently the B(6j) 

bounds will tend to encircle the origin as <a ► oo As 

explained in the earlier section (section 331) this is 
undesirable as it could mean an unrealxsable G(s) 
Consequently, for strictly proper plants, as a necessary 
condition for the realisability of G(s) eCmax) should be 
limited to 3 CH0R76] , CH0R79] 

333 Uncertain unstable plants 

Let "np" be the total number of unstable poles of the 
system P(s) given by equation 22 If each unstable pole 
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remains in the complex right half plane (RHP) over the entire 
range of uncertainty, then this implies that np is a 
constant Consequently, for stability, assuming no pole-zero 
cancellations in P(s)G(s), the fuzzy loop gain polar plot 
should encircle the -1+J0 point np times in the appropriate 
sense The fact that np is a constant removes any ambiguity 
regarding the encirclement count 

However, this is not the case when there are certain 
unstable poles which become stable i e migrate to the left 
half plane (LHP), for certain parameter values In this case 
np and hence the Nyquist encirclement count become variables, 
complicating in the process the stability issue No attempt 
has been made in this thesis to address this problem of a 
variable np A problem with a constant np (see example 5 2, 
chapter 5), however, has been tackled in this thesis. 

As an interesting aside, however , the problem of a 
variable np can be looked upon as a simultaneous 
stabilisation problem of the following kind - 

Let P(s,q) = Pi(s,q) ^ (3 12) 


where Pi(s q) is stable and each ai e [a~ , a^J such that 
a” i 0 or a” < 0 and a^ > 0, i e at could be a 

It t 

non-migratory (i.e will always remain in the RHP) or a 
migratory (i e can cross over to the LHP) unstable pole 
Now define the following two sets of systems — 


P8(s,q) = Pi(s,q) 


1 0 


f| (s - at) 


3 13(a) 


at € Ca , 


a^ < 0, 


where 


0 ). 
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and 

Pus(s,q) = Pi(s,q) > 3 13(b) 

pf ( s " au) 

L 

where at e [0 , a"^] , a*^ > 0 , or 

If L 

ai, € ta f ^ 3 9 a >0 

V V L 

It IS fairly obvious that a 6 ( 5 ) which is admissible 
performance and stability wise, for both P 9 (s,q) and Pu 8 (s,q) 
with the same set of performance specifications as for the 
given system P(s,q), will also be a solution for the original 
problem The point to be noted here is that in this approach 
as Pu 9 (s,q) has a fixed np and Pa(s,q) is anyway stable, the 
issue of stability can be handled relatively easily, as there 
will be no ambiguity regarding the Nyquist encirclement 
count 


The problems associated with the computational aspects 
of the QFT design algorithm described in this chapter and 
consequent remedial measures adopted, constitute the topic for 


discussion in the next chapter 
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CHAPTER 4- 

ALTERNATIVE QFT ALGORITHM 

In this chapter Me Mill discuss some problems associated 
with the QFT algorithm described in section 3 2, chap-ter 3 A new 
algorithm, which takes into account only the extremal open loop 
behaviour and which forms the basis for the QFT design software 
developed in this thesis, is described in this chapter 

4.1 Alternative QFT design algorithm 

The QFT algorithm described in section 3 2, chapter 3, has 
the following disadvantages 

(1) The boundary of the plant template defines the worst 
case open and closed loop behaviour Consequently, an 
accurate generation of the template boundary is 

essential for reliable design Unfortunately, however, 

accurate and computationally efficient algorithms for 
generation of the template boundary do not exist, except 
for some simple cases CBAI883, CBAI89] , CFU90 > , CBART93] 
The most commonly used method is the grid method Here, 
each uncertain parameter interval is partitioned or 
gridded in some suitable manner and at each grid point 
in Q (refer equation 21, chapter 2) the system 
frequency response (magnitude and phase) is computed 
Unless the grid is extremely fine, this method usually 
results in an inaccurate definition of the template 
boundary which can lead to unreliable designs [BART93] 
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(2) Increase in accuracy of the design procedure can be 
achieved only with the computation of a large number of 
values for the B(aj) bounds (refer sec 3 2, chapter 3) 
and also for the templates as well The consequent 

computational labour and memory storage problems could 
be quite phenomenal and will get further compounded with 
increase in the number of uncertain parameters This 

could also possibly force a reduction in the number of 
trial frequencies that can be considered, thereby 
impairing the reliability of the design 

Now a little contemplation will show that the following 
points stand out with regard to the computational aspects of this 
algorithm - 

(a) B(6>) bounds - These bounds are needed only if the 

Nichols chart is used for the design Having a computer 

program to check for satisfaction of the performance 
specifications (equations 2 5, 3 8) should dispense off 
with the computation of these bounds 

(b) Template formation - The actual template, at each trial 
frequency ojo, is bounded b/ the region ABCD(<»)o) , shown 
in the figure below 



Figure 4 1 The ABCD(«o) template 
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This region ABCD(a>o), is obtained by bounding a region 
in the complex plane by the extrema of the loop gain 
magnitude and phase It is obvious that a design for 
the given system, assuming it to have the ABCDCcoo) type 
of uncertainty, will hold for the actual uncertainty as 
well We may thus solve the design problem assuming the 
ABCD(<a)o) type of uncertainty, dispensing off, in the 
process, the generation of the actual templates or, for 
that matter, the entire boundary of each template 

On the basis of these observations, we decided to replace the 
computation of each template (entire or its boundary alone) and 
the B(£o) bounds by the computation of the following values - 

(1) The loop gam extrema (of the magnitude and the phase) 

(2) The worst case extrema of the concerned system transfer 

function (e g Tyr(s), Tydi(s), etc ) to check for 

satisfaction of the performance specifications 

This change yields the following advantages ~ 

(1) Only four values, viz the loop gam extrema, have to 
computed for each template Also, the B(6)) bounds are 
not needed The resulting savings in computational 
labour and memory space could be phenomenal 

(2) Considerable computational savings made, as mentioned 
above, also enable a substantially larger number of 
trial frequencies to be considered, thereby improving 
the reliability of the design 

(3) Explicit computation of the extremal open loop 
behaviour helps m visualising a rough open loop polar 
plot without resort to any graphics utility 


Stability, 
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can then be quickly determined Thus dependence on 

graphical utilities, which is considerable in the 

earlier approach, is greatly reduced here, thereby 
enhancing the user friendly aspect of any algorithm 
based on this approach 

An improved QFT algorithm based on this approach is as 

f ol lows 

STEP 1 Time domain specifications are converted to 
the frequency domain on the basis of a semi - 
heuristic approach described in section 2 2, 
chapter 2 

STEP 2 Define an appropriate, finite set of trial 
frequencies Q 

Assume some appropriate realisable transfer 
function as an initial guess for G{s) 
Preferably, assume G(s) =10, as it will 

give a rough idea about the amount of gam 

needed for satisfaction of the performance 
specifications 

STEP 3 Check if 3(5) is admissible performance wise 
at each tr al frequency If 3t some trial 
frequency <«)o S(s) is inadmissible, then a 
a bound (lower or upper, depending on the 
design problem) on the allowable |S(j<*>o)j, is 
returned Such bounds can be used to define 
another G(s) candidate and this step is then 
repeated 

Check if G(s) is admissible, 


STEP 4 


stability wise 
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If it IS not, assume another structure for 
G(s) and go back to step 3 

STEP 5 If the pre-filter is not needed, then end the 
design procedure, else synthesise a minimum 
phase transfer function F(s) which satisfies 
equation 3 9 and end the design procedure 
The resulting G(s) and, if designed, FIs) will then 
constitute a solution for the QFT design problem This 

approach, for want of a better name, is called as the mini-max 

approach 

A disadvantage, however, of this design procedure could be 
overdesign i e unnecessary design For example, assume a minimum 
phase system with a ABCD(«o) type uncertainty at the real axis 
crossover as shown below 



Figure 4 2 Possible overdesign in the mini-max approach 

Observe that the actual fuzzy open loop polar plot avoids 
encircling the -1+J0 point thereby indicating a stable closed loop 
system However, the ABCD(a)a) template at the real axis crossover 
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indicates an unstable system as it includes the -l+j0 point 
Consequently! based on this erroneous information one might go for 
lead or lag compensation as deemed fit which, however, in 
actuality is unnecessary 

But it may be noted that if the width and/or breadth of the 
ABCDIcm) template at the real axis crossover, crucial for 
stability considerations, is quite small, then the resulting 
overdesign will also be correspondingly less The advantages of 
this approach, mentioned earlier, will then certainly outweigh the 
disadvantage of overdesign Consequently, we may, under such 
conditions, overlook such overdesign 

4.2 Computational problems associated with the mini-max approach 

Ex tremization of the loop gain and the concerned system 
transfer function (e g Tyr(s), Tyd2(s), etc ) is a very crucial 
step in the mini— max approach This problem of ex tremisation , 
a non-linear programming problem CRA0B4] , can be stated 
mathematically as - 

Extremise f(q), q e Q (41) 

where f(q) could be the magnitude or pha'"e of the specified 
transfer function and q is the uncertain parameter vector (refer 
equation 2 1) 

However, the solution of equation 4 1, obtained using any 
efficient numerical nonlinear programming method is always a local 
extremum which need not be global over Q CRA084]. Consequently, 
there is the danger of ending up with erronous results and hence 
erronous designs 


This fact forced us to take 


up the following 


two 
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computational alternatives - 

(1) One alternative is to carry out a blind search over the 
entire domain of Q for the extrema Computationally, it 
means that we arbitrarily partition Q into a grid gr(Q) 
and at each grid point q* e gr(Q), compute f ( q* ) A 

straight forward search among the resulting f(q*) values 
will then give an approximation to the actual extrema 
The accuracy of this method, obviously, will depend on 
the fineness of the grid gr(Q) Increase in accuracy 
however, can be achieved only with an increase in 
computational labour This method, known as the grid 
method, also becomes computationally inefficient and/or 
intensive with increase in the number of uncertain 
parameters This method, however, at the worst, is as 
inefficient as the earlier given algorithm (section 3 3, 
chapter 3) 

A mini-max algorithm based on the grid method with 
uniform partitioning of the grid has been implemented in 
this thesis 

(?) Another alternative is based on an a prion knowledge of 
a small region R in Q, where the extrema are guaranteed 
to be located The region R should be small enough so as 
to search and locate the extrema with the minimum 
computational labour and maximal accuracy possible 
Powerful optimisation routines like those found in the 
NAG or IMSL math libraries can be utilised for this 
purpose This method, if possible to apply, would be 
more efficient and faster than the grid method, as shown 
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by the comparison between these two methods in sec. 4.3,2. 
As the extrema search is carried out in a restricted 
region R, this method can be called as the restricted 
region method . 

However, no analytic means exist for finding this 
region R, except in the case of a class of uncertain 
systems known as linear affine uncertain systems. The 
restricted region method, as applied to this class is 
described in the following section. 

4.3 Restricted region method as applied to the class of linear 
affine uncertain systems 

In this section we will be needing the following definitions 

Definition 4.3.1- A function farlR" »• K, is said to be 

linearly affine if it can be expressed as a 

sum of a linear function grR" ► R and a 

constant k e R, i.e. 
fa(q) = g(q) + k , q s K”* 

Definition 4.3.2- Any system with a structured uncertainty 

model P{s) giyen in equation 2.2, 
chapter 2, is said to be a linear affine 
uncertain system, if each model coefficient 
ai and bv, are linearly affine in q, the 
uncertain parameter vector. 

Definition 4.3.3- Given the n - dimensional rectangle 

Q = £q e R"lq~<q <q'^;l<i^n} 
an edge of Q is defined as the set of 


vectors Ei S Q, given as 
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EL = {e € Q |e~<e <e* and 

' V U (. 

e = e~ or e'^ , j ^ i } . 

J J J 

Now, as G(s) and F(s) are rational transfer functions, it is 
very easy to prove the following lemma — 

Lemma 4.3.1 - CBARri93] If the uncertain system P(s) shown in 
figure 2.1 is a linear affine uncertain system, 
then all the system transfer functions for this 
feedback configuration (e.g. Tyr(s), the 
complementary sensitivity function, etc.) are 
also linear affine uncertain systems. 

Consequently, henceforth, we will concern ourselves with some 
general linear affine uncertain system, viz. V(s) with the 
knowledge that whatever results hold for V(s) will also hold for 
all the system transfer functions (including the open loop gain) . 

Now use will be made of the following very important theorem 
proved in CFU90] . 

Theorem 4.3.1 - For a linear affine uncertain system V(s), at 

each trial frequency coo, the V(jcoo,q) 
template boundary is the image of all or some 
edges Ev of J. 

The implication of this theorem is that, for linear affine 
uncertain systems, the frequency response extrema at each trial 
frequency lie on the edges of Q. The set of edges of Q, will then 
constitute the restricted region R, mentioned in the earlier 

section . 

Extrema search along each edge Ei of Q, is essentially a 
one dimensional search as is evident from definition 4.3.3. There 
are very many one dimensional search algorithms CRA084], of 
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which we chose the quasi Newton optimisation technique. Greater 
speed and accuracy, and the fact that it is the easiest to 
implement as far as NAG optimisation routines CNAG90] are 
concerned, were the reasons for choosing this particular 
optimisation technique- More specifically, a mini-max algorithm 
using the NAG routine E04JAF, an easy to use quasi Newton 
algorithm CNAG90], was developed in this thesis for handling 
linear affine uncertain systems. The extrema search is restricted 
to the edges of Q. 

4.3.1 A Digression 

A fairly recent paper by A. C. Bartlett CBART93] , however , 
has shown that an extrema search of the kind just described 
above, for linear affine uncertain systems, is really 
not needed, as far as magnitude ex tremisation is concerned. 
Here a very simple method for finding the magnitude extrema 
using at most four values of V(j<oo,q) along each edge Et of 
Q, IS given. We, however, persisted with the use of the NAG 
based extrema search algorithm for the following reason - 

We observed that our algor thm,on an average ,made about 
20 f(q) value computations ts against 5 (worst case) 
made by Bartlett s method. In other words, our 
algorithm, on an average, executed 5 times more the 
worst case number of arithmetic operations as executed 
in Bartlett s method. Now, even in a PC an increase by 
a factor of this magnitude in the number of arithmetic 
operations would be insignificant, computation time 
wise. So, efficiency wise, our algorithm can be said to 
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more or less as efficient as Bartlett s method. 


4.3.2 Performance comparison of the restricted region method 
and the grid method 

The test system chosen for this comparison was the one 
handled in example 5.3, chapter 5, viz. 


P(s) = 


(s + a) 


s® + (13 + b)s^ + (30 + 13b)s + 30b 


where a e [0.5, 2. 5] and b e [4. 0,8.0]. 

For the grid method each uncertain parameter s interval was 
uniformly partitioned into 99 sub-intervals. 

Frequency response computation for the system sensitivity 

1.0 


S(s) = 


1.0 + L(s) 


over a set of 100 frequencies belonging to the interval 
[0.01,100] rad/sec. resulted in the following observations - 



Grid 

Restricted Region 


method 

method 

Time taken 
for each 
frequency 

8 

minutes 

12 

seconds 

(on an average) 




Table 4.1. Comparison of the performance of the grid and the 
restricted region methods. 

As expected the restricted region method was found to 
be quite fast relative to the grid method. These 
computations were done on a HP 9000 (series 800) super 
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mini computer. 

4.3.3 Beyond the linear affine uncertainty case 

Theorem 4.3.1 is the most advanced result concerning 
frequency response computation for systems with 
structured uncertainty CBART933. For uncertain systems 
with functional dependence of the coefficients more 
complex and nonlinear than the linear affine case, to 
date, no analytical result exists which might ease 
the computational aspects. It is interesting to note 
the existence of a similar situation in the area of 
Kharitonov - like analysis of parametric uncertainty 
systems CBARM933. 

In view of the existing scenario, for uncertain 
systems with uncertainty more complex than that 
encountered in the linear affine case, the grid method 
(section 4.2) seems to be the only alternative. 

4. 4 Nature of structured uncertainty handled by the software 
developed in thi*^ thesis 

We will need the following definition in this section. 

Definition 4.4.1 - An uncertain system with a structured 

uncertainty model given in equation 2.2, 
is said to be polynomial ly uncertain 
if each coefficient (a». or bv) of the 
model in equation 2.2 is some 
multivariable polynomial in q, the 


uncertain parameter vector, with real 
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coefficients . 

Note that a linear affine uncertain system is a special case 
of polynomial ly uncertain systems. 

The software developed in this thesis is equipped to handle 
arbitrary polynomial ly uncertain systems. The motivation for 
adoption of capability to handle this type of uncertain systems 
stems mainly from the following reasons - 

(1) Most of the systems encountered e.g electrical, 
mechanical, etc., have models with coefficients which 
are some multivariable polynomials of the system 
parameters CTHA60]. 

(2) For each coefficient, say av, such that 

av = fn(q) , q € Q 

where fn( ) is some non - polynomial, nonlinear 

continuous function of q, it is always possible, by the 
Weierstrass approximation theorem [SIMiSl], to find a 
polynomial g(q) with real coefficients such that 
I f(q) - g(q) I < £ f over Q 

where £ > iS , is an arbitrarily small number. If £ is 
sufficiently small then we can take 
ai = g(q) , over Q 

and carry out the design for the resulting "modified" 
system. If the approximation effected is quite accurate 
then this design should hold for the actual system as 
wel 1 . 

It may be noted that variable time delay systems can also be 
handled by the software developed, by going for rational Pad6 
approximations for the " e " terms involved. 
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The solutions to five examples solved using the algorithm 
described in this chapter (section 4.1) are presented in the next 
chapter. 
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CHAPTER 5 

EXAMPLES 

In this chapter we will discuss the solutions to five design 
problems, obtained using the mini-max approach to QFT described in 
the last chapter (see section 4.1, chapter 4). 

Here, we will be needing the following definition - 
Definition 5.1 - A frequency interval W = [oiijoiz] is said to 

lagspaced into n sub-intervals if it is 
partitioned into n sub-intervals each of 

length Clog ( < 01 / <<> 2 ) ] / n . 

10 

NOTE : This feature was adopted from a MATLAB 
utility called "logspace" CMATS73. 



where "col" is the motor shaft speed and "e" -■5 the motor 

armature voltage and where 

ai = 0.056 X 10'^ + 0.4(J1) 
az = 3 X 10"'* + 22 X 10“'‘(J1) 

The load moment of inertia J1 is assumed to vary in the 
interval [0.7x10 f 14x10 kg— m*. 

The performance goal of interest here is to have the 
motor shaft speed track a speed reference signal "r". 
Response speed and damping considerations lead to the 
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Tol lowing time domain bounds to be specified - 

SPl - 2'/. settling time of the unit step response 

ts € C6.0,a.0] seconds. 

SP2 - Peak overshoot of the unit step response Mp < 207. . 

NOTE : The assumption of a lower bound on ts allows us to set up 
the problem, relatively easily, as one in which the frequency 
domain specifications are of the form in equation 2.5(a). Even 
though this specification may seem unrealistic, it should be noted 
that our main objective here was to illustrate the use of this 
software rather than solve a real world design problem. 

Solution - The design procedure consisted of the following 


steps - 

STEP 1 Setting up of the frequency domain specifications - 

The time to frequency domain conversion needed for this 

problem is described in example 2.1, chapter 2. The 

resulting bounding transfer functions, as found there, are 

(0.216s + 1.08) 


Tyr(s) = 


s + 0.94s + 1.08 


yr 


(S) = 


1.35 


+ 1.998s + 2.736s + 1.35 


The equivalent performance specifications in the frequency 
domain are now to design a G(s) and a F(s) to ensure 

I T yr (jco)| ^ | Tyr(ja))| S | Tyr ( J to ) | , V<o ^ 0 (5.1) 

(see section 2.2, chapter 2). 

STEP 2 Choice of the trial frequencies set Q and an initial 
guess for G(s) - 

O was chosen to be a set of 100 logspaced frequencies 
(see definition 5.1) in the (arbitrarily chosen) range of 
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[0.1,1000] rad/sec. 

G(s) = 1.0 Mas the initial guess for G(s). 

STEP 5 Synthesis of G(s) — 

The initial guess of G(s) = 1.0 was found to be 
satisfactory in that it ensured 

ALm[Tyr( J0>) ]Tn<ix < Lm[Tyr(jO))] — LmCTyr( J<*>) ] , 

Voj € O (5.2) 

It also resulted in a stable closed loop set-up. This 
particular fact can be verified using the Routh-Hurwitz test. 

Thus simple closing of the loop around P(s) was found 
to be satisfactory, performance and stability wise. 

STEP 4 Synthesis of the pre— filter F(s) - 

Bounds on the magnitude of the pre-filter F(s) were 
computed at all the trial frequencies, using equation 3.9, 
chapter 3. Use of these bounds and MATLAB lead to a design 
of 

F(s) = , for the pre-filter. 

s* + 1.3s + 1.0 


val idated 
figure 5. 
parameter 

by J1 IS 


Time domain testing of this design using MATLAB 
it as can be seei from the unit step responses ,in 
1, obtained for a s ’■t of 10 values of the uncertain 
J1 (also see table 5.1). 

The satisfaction of equation 5.1, for all values taken 
shown in figure 5.2. 

Stability can be established via the fuzzy open loop 


polar plot, shown by means of a finite set of discrete ABCD(a>o) 
type templates (see section 4.1, chapter 4) in figure 5.3. 



Parameter 

value 

Settl ing 

time (sec) 

Peak 

overshoot {’/.) 

0 . 7e-4 

6.03 

7 

i.4e-4 

6.03 

7 

2.8e-4 

6.03 

7 

35e-4 

6.03 

7 

54e— 4 

6.03 

7 

77e-4 

6.03 

7 

91e-4 

6.03 

7 

10e-3 

6.03 

7 

12e-3 

6.03 

7 

14e-3 

6.06 

7 


Table 5.1 Settling time and peak overshoot of the 
unit step response for different values 
of the uncertain parameter J1 - Example 1. 


5.2 Example 2 

Problem Statement - The output of an unstable uncertain 
system 


P(s) = 


a « [1.0, 5.0] 


' (s - a)(s + 10) ’ - 

has to satisfy tracking specifications same as those given 

in example 5.1, viz. SPl and SP2. 

Solution - The design procedure consisted of the following 
iteps - 

STEP 1 Setting up of the frequency domain specifications - 
Since the time domain specifications are the same as in 
example 5.1, the equivalent frequency domain specification 
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15, as in example 5.1, to design a G(s) and F(s) so as to 
ensure that Tyr(s) satisfies equation 5.1 with the bounding 
transfer functions Tyr(s) and Tyr(s) being the same as in 
example 5.1. 


STEP 2 Choice of the trial frequencies set O and an 
initial guess for G( 5 ) — 

O was chosen to be a set of 100 logspaced frequencies 
in the (arbitrarily chosen) range [0.1,1000] rad/sec. . 

G(s) = 1.0 was the initial guess for G(s). 

STEP 3 Synthesis of G(s) - 

The initial guess of G(s) = 1.0 indicated that the 
maximum lower bound on | G(jw) |, Vto e Q, was 3565 . 

A candidate G(s) = 3565.0 was consequently found to be 
satisfactory performance wise (i.e. it satisfied equation 
5.2) and also ensured a stable closed loop system. In order 
to, however, provide a "realistic" design involving a 
strictly proper G(s), we added a simple pole to this 
candidate G(s) giving 

G(s) = — ^ , b > 0. 

However to ensure stability we had to take a very large 
value for b (greater than 1000). In order to keep b not too 
large and at the same time ensure stability, we additionally 
resorted to a lead compensator resulting, for b = 300, in the 

following design for G(5) 

34.224 X 10®(s + 25) 

(s+90)(s+ 300) 


This 6(s) was found to satisfy equation 5.2 as well as 
ensured closed loop stability. 

STEP 4 Synthesis of the pre-filter F(s) - 
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Time domain testing of this design using MATLAB 
validated it, as can be seen from the unit step responses in 
figure 5.4, obtained for a set of ten values of the uncertain 
parameter "a" (also see table 5.2). 

The satisfaction, by this design, of the frequency 
domain specification of equation 5.1 is shown in figure 5.5. 

Stability can be established from the fuzzy open loop 
polar plot, consisting of a finite number of discrete ABCD(<oo) - 
type templates and shown in figures 5.6 and 5.7. 


Parameter 

values 

Settling 

time ( sec . ) 

Peak 

overshoot ( ) 

1.0 

5.34 

9.7 

5.0 

5.34 

9.44 

1.5 

5.34 

9.54 

4.5 

5.34 

9.6 

2.1 

5.34 

9.35 

4.1 

5.34 

9.72 

3.2 

5.34 

9.01 

3.5 

5.34 

8.92 

1.25 

5.34 

9.61 

4.7 

5.34 

9.54 


Table 5.2 Settling time and peak overshoot of the 
unit step response for different values 
of the uncertain parameter a -Example 2. 
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5*3 Examplo 3 

Problem Statement - Given the system 

P(S) = (s 

(s+3)(s+10)(s+b) ’ 

a € [0.5, 2. 5], b € [4. 0,8.0] 

find a stabilising compensator G(s) such that the sensitivity 

' 1 0 

function S(s) = satisfies 

1.0 + L(s) 

I S(jco) I < 0.01 , V CO € [0,0.3] rad/sec. 

Solution - The design procedure consisted of the following 
steps - 

STEP 1 Setting up of the frequency domain specifications - 
The design objective was to design a G(s) to ensure 

I SCjoj) I < 0.01 , Voj e [0,0.3] - given 

and I S(joj) I < 100 , V<o > 0.3 - this 

specification was arbitrarily 
made keeping in mind that 
I S ( J CO ) I mooc >1.0 (5.3) 

STEP 2 Choice of the trial frequencies set Q and an 
initial guess for GCs) - 

O was chosen to consist of a set of 100 logspaced 
frequencies in the (arbitrarily chosen) range of [0.01,100] 
rad/sec. . 

G(s) = 1.0 was an initial guess for G(s). 

STEP 3 Synthesis of G(s) - 

The initial guess of G(s) = 1.0 indicated that the 

maximum lower bound on | G(j<o) | for all the trial frequencies 

was 4320 . Consequently we chose a candidate GCs) 
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G(s) = 


4320(b) 


(s + b) 

"Tuning" the value of b we finally 
(arbitrary choice) resulting in 


b > 0 

?ttled for b 


= 350.0 


G(s) = 


15.12 X 10' 
(s + 350) 


which satisfied the performance specification of equation 5.3 
as well ensured a stable closed loop system- This can be 
clearly seen from figures 5.8 - 5.11. 


5.4 Example 4 

Problem Statement - For the given system 

= (s ^ 2?(s ^ b) - ’ ^ Cl. 0,10.0] 

beset with impulse-like disturbances (like e.g. voltage 

spikes) at the plant input, a compensator G(s) is sought such 
that the system response to these disturbances has a 27. 
settling time of less than a second, for all values of the 
uncertain parameters a and b. 

Solution - The design procedure consisted of the following 


steps - 

STEP 1 Setting up of the frequency domain specifications - 
We assumed unit impulse disturbances at the plant input. 
Consequently, the design objective was to ensure that the. 


impulse response of the transfer function 

y(5) _ P(s ) 


T ydl ( s ) = 


dl ( s ) 


1.0 + L(s) 


(5.4) 


(see figure 2.1, chapter 2) has a 27. settling time of less 
than a second for all values of the uncertain parameters a 


and b . 

To have such a settling time, as a rule of thumb, we 
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define a bounding transfer function Tydi(s) with a settling 
time of less than 1 second, such that ensuring 

I Tydi(ja>) j max < j Tydi(jco) j, Vco S: 0 (5.5) 

would ensure the disturbance response settling time 
specification. This frequency domain specification of 

equation 5.5 was made on the same lines as suggested in 
[D'AZSB] and CJAY92]. 

Consequently we considered a boundary model of the form 


T ydl ( S ) 


2 

con 


2 2 
S + (2^ con) 5 + con 


(5.6) 


It can be easily verified that the 2/C settling time for 

this model is given as 

^ loge(50con) 

ts = = 

< con 

We arbitrarily assumed con = 14 rad/sec (which is far 
greater than the average system bandwidth of about less than 
2 rad/sec.). This fixed up ? > 0.47 for ts < 1 second. As a 
result, we assumed ? = 0.5. Thus the boundary model used was 


Ty<li(s) = 


196 


s* + 14s + 196 


L.-'* > L L'iRaH 

! 1 r KANPUR 


IflL No. 


A.WM 


With a settling time of 0.94 sec. . 

STEP 2 Choice of the trial frequencies set Cl and an 

initial guess for G(s) 

a was chosen to be a set of 100 logspaced frequencies 


in the (arbitrarily chosen) range [0.01,100] rad/sec 
B(s) = 1.0 was an initial choice for G(s). 

STEP 3 Synthesis of G(s) - 


The initial guess of G(s) 


= 1.0 indicated that a gain 
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of more than 1.08 was required over a certain low frequency 
range, essentially at around 3 rad/sec. Based on this a 
controller candidate 


G(s) = 


45.0 


(s + 30) 

was considered which was found to to satisfy equation 5.5 as 
well as ensured a stable closed loop system. 

Time domain testing using liATLAB, however, showed that 
for low values of the parameters a and b (e.g. a = b = 1.0) 
the resulting settling time was, contrary to expectations, 
greater than 1 second . 

Investigation using the root locus indicated that this 
could possibly be due to the relative closeness of the closed 
loop poles to the jto - axis, for low values of the uncertain 
parameters a and b. 

Now, note that the two poles of PCs) can be weaned away 
from the jco - axis by getting G(s) to have two LHP zeros. 
These two zeros were chosen to be located at - 8 + j3. The 
Idea behind this choice was that if the resulting loop 
function L(s) has a dominant pole pair at these locations 
th >n the settling time t* would be less than 1 second. This 
emerges from the fact that a second order system with the 
same model as in equation 5.6 and with the poles located at 
these locations would have a settling time of less than a 


second (specifically about 0.76 sec). 

To make this G(s) realisable we added two "far-off 

poles at - 41 + jl0, resulting in the following design 

292. 8(s^ -*• 16s + 73) 

(s* + 82s + 1781) 


G(s) 
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The gain factor of 292.8 was required to "push" the 
open loop poles towards the zeros at - 8 + j3. 

This G(s) was found to ensure satisfaction of equation 
5.5, as well as closed loop stability. 

Time domain testing for a set of ten value pairs of a 
and b (see table 5.3) also validated this design. The 'resulting 
impulse responses are shown in figure 5.12. 

Satisfaction of the frequency domain specification of 
equation 5.5 is shown in figure 5.13. 

Figures 5.14 and 5.15 show the resulting fuzzy open loop 
polar plot consisting of a finite set of discrete ABCDtoio) - type 
templates which were used in establishing closed loop stability. 


Parameter 
value - a 

Parameter 
value - b 

Settling 
time ( sec ) 

1.0 

1.0 

0.95 

10.0 

10.0 

0.61 

2.0 

9.0 

0.71 

9.0 

2.0 

0.56 

3.5 

6. 5 

0.68 

6.5 

3. 5 

0.6 

4.7 

7.3 

0.65 

7.3 

4.7 

0.6 

5.5 

8.2 

0.64 

8.2 

5.5 

0.6 


Table 5.3 Settling time of the unit impulse response 
for different values of the uncertain 
parameters a and b - Example 4. 
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5.5 Exan^ale 5 

Problem Statement — For the given system 
P(s) = 


(S + 3)(s + 10)(s + b) 
a e C0.5,2.5] and b « C4. 0,8.0] 

beset With impulse-like disturbances at the plant output, a 

compensator G(s) is sought such that the system response to 

these disturbances has a 27. settling time of less than a 

second for all values of the uncertain parameters a and b. 

Solution - The design procedure consisted of the following 


steps 


STEP 1 


Setting up of the frequency domain specifications - 


We assumed unit impulse disturbances at the plant 

output. Consequently, the design objective was to ensure 

that the impulse response of the following transfer function 

y(s) _ 1.0 


T ycl2 ( S ) = 


(5.7) 


d2(s) 1.0 + L(s) 

(see figure 2.1, chapter 2) has a settling time of less than 

a second for all the values of the uncertain parameters a and 


b. 

As in example 5.4, we needed a boundary transfer 
function Tyd 2 (s) with a settling time of less than 1 second 
such that ensuring of the following frequency doi’.ain 
specification 

I Ty<l2( J Co) I TnOLX i [ Tyd2(j0)) | , V<i> i 0 (5.8) 

would, as a rule of thumb, ensure the disturbance response 
settling time specification. 

We, on similar lines as in CD AZ8S3 chose the following 
boundary transfer function model 
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T ydZ ( S ) = 


(s + a) 


■+ a) ^ + f}^ 


(T 


(5.9) 


a 


with the 2X settling time given as 


ta - + /? )/a] 

ot 


For a settling time of less than 1 second we chose 
( arbitrari 1 y ) a = 10.0 and ft = 50.0, resulting in the 

following boundary transfer function 


TydzCs) 


260. 0(s 10) 

s^ + 20s + 2600 


with a settling time of 0.95 sec. 

STEP 2 Choice of the trial frequencies set Cl and an 
initial guess for G(s) - 

O was chosen to be a set of 100 logspaced frequencies in 
the (arbitrarily chosen) range of [0.01,100] rad/sec. 

G( 5 ) = 1.0 was the initial choice for G(s). 

STEP 3 Synthesis of G(s) - 

The initial guess of G(s) = 1.0 was found to be 
satisfying equation 5.8 and also ensured a stable closed 
loop system. But MATLAB aided time domain testing of this 
design sho 'ted that for low values of a and b, the resulting 
settling time was, contrary to expectations, greater than 1 


second . 


Investigation using the root 

locus, as 

done 

in 

example 

5.4, indicated that 

this could 

possibly 

be 

due 

to the 

relative closeness of 

some closed 

loop poles 

to 

the 

jco - axis 


for low values of the uncertain parameters a and b. 

Based on this investigation, another candidate G(s) 
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given as 

G(s) = 9-58(s + 3.5) 

( s + 30 ) 

was c on s 1 d sr"sd . The zeno at ~3.5 was added to "take cane" of 
of the open loop pole at -3.0. Note that due to the 
non-minimumphase open loop zero at a, we have the 
complementary root locus CKU091]. The pole at -30.0, of 6(s) 
was simply a "far-off" pole added to make 6(s) realisable. 
Observe that this 6(s) is a lead compensator. 

This G(s) was found to satisfy equation 5.8 and also 
ensure a stable closed loop system. This design was also 
validated by MATLAB aided time domain testing for a set of 
ten value pairs of the uncertain parameters a and b (see 
table 5.4). The resulting impulse responses are shown in 
figure 5.16. 

The satisfaction of the frequency domain specification 
of equation 5.8 is shown in figure 5-17. 

The resulting fuzzy open loop polar plot consisting of a 
finite set of discrete ABCD(ciio) - type templates is shown in 
figure 5.18. 


Parameter 
value - a 

Parameter 
value - b 

Settling 
time ( ~ec ) 

0.5 

4.0 

0. SS 

2.5 

8.0 

0.52 

ID 

• 

8.0 

0.47 

2.5 

4.0 

0.77 

1.5 

6*0 

0.58 


Table 5.4 


continued on the next page 
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Parameter 
value - a 

Parameter 
value - b 

Settl ing 
time (sec) 

1.0 

7.0 

0.52 

2.0 

5.0 

0.66 

1.8 

4.5 

0.69 

2.2 

6.2 

0.59 

0.8 

7.8 

0.48 


Table 5.4 (cantd.) Settling time of the unit 
impulse response for different 
values of the uncertain parameters 
a and b - Example 5. 


Thus in this chapter we described the solutions to five 
examples, obtained using the mini-max approach based QFT design 
scheme described in the previous chapter . 

Though the semi-heuristical conversion of time domain 
specifications to the frequency domain worked in the first two 
examples, it however ran into trouble in the last two examples. 
One reason for this could have been improper choice of the 
boundary transfer functions. This particular aspect, however, 
uTfortunately could not be investigated further. But t ->e last two 
ex 'mples certainly sound a cautionary note regarding thw setting 
up of frequency domain specifications equivalent to a given set of 
time domain specifications. 


r 
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FIGURE 5.4 STEP RESPONSE EXAMPLE 2 
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FIGURE 5.5 BODE MAGNITUDE PLOT OF ALLOWED AND OBTAINED CLOSED LOOP GAIN - EXAMPLE 2 
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IMAG 



FIGURE 5.7 POLAR PLOT (MAGNIFIED) OF OPEN LOOP GAIN EXAMPLE 2 
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FIGURE 5.11 STEP RESPONSE EXAMPLE 3 
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MAGNITUDE (dB) 



FIG. 5.13 BODE MAGNITUDE PLOT OF ALLOWED AND OBTAINED DISTURBANCE - EXAMPLE 4 



REAL 


IMAG 
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frequency polar plot of open loop gain EXAMPLE 4 




REAL 



FIGURE 5.15 HIGH FREQUENCY POLAR PLOT OF OPEN LOOP GAIN -• EXAMPLE 4 
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FIGURE 5.16 IMPULSE RESPONSE EXAMPLE 5 
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FIG. 5.17 BODE MAGNITUDE PLOT OF ALLOWED AND OBTAINED DISTURBANCE - EXAMPLE 5 


REAL 
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IMAG 



FIGURE 5.10 POLAR PLOT OF OPEN LOOP GAIN EXAMPLE 5 
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CHAPTER 6 

CONCLUSIONS 

An issue in feedback design, which has been a subject of 
intense research in the last few years, is that of ensuring system 
performance in the face of large variations in the values of the 
system parameters. That QFT is quite effective and easy to apply 
in solving this design problem, at least for the SISO case, has 
been demonstrated in this thesis through the solution of five 
varied examples (refer chapter 5). A software package, quite 
general enough to handle arbitrary polynomial ly uncertain systems 
(refer section 4.4, chapter 4), has been developed to provide 
computer — aided QFT design. This software has the following unique 
features - 

(1) The use of the Nichols chart, as advocated in classical 
QFT CJAY93], is dispensed off with. The resulting 
saving in computational labour is phenomenal. 

(2) Dependence on graphical utilities, which is very high 
in classical QFT [BAI91], CBAR93] , is greatly reduced, 
when using this software (refer appendix A). It also 
reveals QFT to be essentially a natural extension of the 
traditional loop shaping based design of the kind found 
in e.g.CBOWSB]. These factors help in enhancing the user 
friendly aspect of QFT. Anyone, with a basic knowledge 
of classical loop shaping, can quickly and effectively 
learn to grapple, using this software, with realistic 
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design problems involving large parametric uncertainty. 

6.1 Beyond single loop LTI SISO QFT design 

Though only continuous time, LTI SISO systems in a single 
loop feedback configuration' (refer fig. 2.1), have been considered 
in this thesis, it may well be noted that QFT has been used to 
handle the following types of systems as well - 

(1) Multiple loop continuous time, LTI SISO systems. 

(2) Discrete time LTI SISO systems. 

(3) MIMO LTI systems. 

(4) Time varying systems. 

(5) Nonlinear systems. 

An exhaustive survey of the literature dealing with these 
design problems is given in CH0R91]. Related references which 
have come out fairly recently and which provide for a better and 
simplified elucidation of the QFT design procedure are - CYAN93] , 
for discrete time SISO LTI systems and CYAN923* CPER93] , for MIMO 
LTI systems. 

6.2 Scope for future work 

6.2.1 Defining performance specifications 

QFT, as described in section 2.2, chapter 2, requires 
time domain specifications to be converted to the frequency 
domain. Lack of a theoretical framework under which this 
conversion can be rigourously effected, forces us to adopt 
semi-heuristic techniques of the type exemplified in example 
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2.1, chapter 2. That this could be problematic is shown in 
examples 5.4 and 5.5, chapter 5. Consequently, this aspect 
of time domain to frequency domain conversion is one area 
where research is definitely required. The theory of 
wavelets CCHI92], needs looking up in this particular regard. 

Fairly recently, however, there have been efforts to 
ensure tracking specifications of the form given in equation 
2.4, directly in the time domain CBAR93]. This approach, 
known as time domain QFT , however, is in its infancy and 
provides another area for further research. 

6.2.2 Usability of the software developed in this thesis 

The use of two NAG routines (refer appendix A) by the 
software developed in this thesis, limits its use to computer 
systems loaded with the NAG math library. Also, a convenient 
and transparent interface to NATLAB has not been provided 
here. Thus there is scope for improvement of this software 
to increase its portability and ease of use. 
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CASH82] 


CAZV92] 


CBAI88] - 


CBAI89] - 


CBAI91] 


CBAR933 - 


CBARM93] 


CBART933 


CBHA873 


CBOD453 
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APPENDIX A 

Here, a brief description of the software developed in this 
thesis IS provided with the help of an example. 

A. 1 General description of the OFT design software 
A. 1.1 Uncertain system specification 

The uncertain system should be specified in the form 
given in equation 2.2, chapter 2. Here, each coefficient ai.(q) is 
assumed to be some multivariable polynomial in the uncertain 
parameter vector q, with real coefficients, i.e.. 


a I. ( q) = Z a. t 
j=o J J 


at 


o 


R, 


n- 1 S 

where t^*’ = n (q ) ^ , 

J k=o k 


where "n" is the number of uncertain real system parameters 
q € Cq~ , q*] S R, while "r" and "s '* are some non-negative 

k k k jk 


integers . 


Each "t^*'" IS known as a "term" in the coefficient at(q) 


A. 1.2 Specification of the trial frequencies' set Q 

The trial frequencies' set O can be specified to be any 
arbitrary finite set of real frequencies, in the units of radians 
per second. Provision is also made to specify a set of logspaced 
frequencies (see definition 5.1, chapter 5) in the interval 
[10'^‘, where di and da are integers with di < da. 


9 



84 


A. 1.3 Problem type specif ication 

The following types of design problems can be handled, 
one at a time, by this software - 

TYPE 1 The performance specification of equation 2. 5Ca3, 
chapter 2, is handled here and ^involves synthesising a 
G(s) and, if necessary, an F(s). Examples 5.1 and 5.2, 
chapter 5, are of this problem type. The performance 
admissibility of | G{jo>) | at each trial frequency ci> e £1 
(see section 3.2, chapter 3) is tested. If inadmissible 
a bound for | G(jto) | is found. This bound will be a 


lower one for a 

non-zero 

Tyr(s) and 

upper 

one for 

Tyr(s) =0. A 

Tyr ( S ) = 

0 allows 

one to 

have a 


performance specification of the form 
I Tyr(j<o)| < </(«"') , V<o > 

where <oh is some high frequency value and n is some 
integer with K > 0. 

TYPE 2 A disturbance attenuation specification of equation 
2. 5Cb>, chapter 2 (see figure 2.1, chapter 2, and 
example 5.4, chapter 5) falls under this design problem 
type. It <nvolves synthesis of G(s) alone, as the 
pre-filter K(s) is not needed here. A lower bound for 
I G(j«) I IS found, if G(s) is found inadmissible at the 
trial frequency e £1. 

TYPE 3 A disturbance attenuation specification of equation 
2. 5Cc3, chapter 2 (see figure 2.1, chapter 2 and example 
5.5, chapter 5) or a sensitivity specification of the 
kind given in example 5.3, chapter 5, falls under this 

It involves, as in a type 2 


design problem type. 
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problem, synthesis of 6(s) alone, as the pre-filter is 
not needed. A lower bound for | G(j(o) | is found, if 
G(s) IS found inadmissible at the trial frequency w e n. 

A.i.4 Testing closed loop stability 

For a performance admissible G(s), at each trial 
frequency oj « O, the loop gain extrema (i.e. of the magnitude and 
phase of G(j&j) ) are found. For a linear affine uncertain system 
(definition 4.3.2, section 4.3, chapter 4), the restricted region 
method (section 4.2, chapter 4) can be used. For any other 
arbitrary polynomially uncertain system (definition 4.4.1, section 
4.4, chapter 4) use the grid method (section 4.2, chapter 4). 
These extrema can be then used to visualise a rough open loop 
polar plot involving the ABCD(ojo) - type templates (section 4.1, 
chapter 4, and e.g. of the form given in figure 5.3, chapter 5). 
Applying the Nyquist criteria to this fuzzy open loop polar plot, 
stability can then be determined without resort to the use of any 
computer graphics utility. 

A. 1.5 Operational details of the software developed 

The programmii g has been done in C and Fortran, with a 
transparent interface between the C and Fortran modules. The 
Fortran module was necessitated due to the use of two NAG routines 
VIZ. CNAG903 “ 

(1) E04JAF - a quasi Newton optimisation algorithm, on 
which the restricted region method is based. 

(2) C02AFF — an algorithm to find all the roots of a 
complex polynomial, which aids phase computation. 
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Actual operational details of the software developed, 
are provided in a help file called "QFT_help". 


A. 2 Exanqile 

Consider the design problem given in example 5.3, chapter 5. 
This IS a type-3 problem (section A. 1.3) with the following 
details, which were input, aided by prompting of the computer 
program - 


Uncertain system - 


PCs) 


(s + a) 

(s + 3) (s + i0)(s + b) 


- (s a) 

s® + (13 + b)s* + (30 + 13b)s + 30b 


Bounds on the - 
uncertain parameters 


a e C0.5,2.53 
b e [4.0,S.0] 


Problem type - This is a type 3 problem with the upper bound 

TydaCjto) =0.1 , V« € [0.0,0.33, and 

Tyd2(j<0) = 100 , Vco > 0.3 

Trial frequencies - O was chosen to be a set of 100 logspaced 

“*2 2 

set frequencies in the range [10 , 10 3 rad/sec. 


Candidate G(s) 


- A candidate G(s) was specified to be 

1 s. 1 “7 V 10® 

6(s) = 


(s 


350) 


This candidate G(s) eventually turned out to be 
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a solution for the given design problem. To see 
the ensuring of closed loop stability by this G(s), 
the loop gain extrema are provided in the following 
table, VIZ. table A.l. 


FREQUENCY 

( rad/sec ) 

0.010000 
0.010975 
0.012045 
0.013219 
0.014508 
0.015923 
0.017475 
0.019179 
0.021049 
0.023101 
0.025354 
0.027826 
0.030539 
0.033516 
0.036784 
0 .040>i70 
0 . 0443L-6 
0.048626 
0.053367 
0.058570 
0.064281 
0.070548 
0.077426 
0.084975 
0.093260 


MINIMUM 

LOOP 

MAGNITUDE 

(dB) 

19.086528 

19.086871 

19.087284 

19.087781 

19.088380 

19.089102 

19.089971 

19.091017 

19.092277 

19.093794 

19.095621 

19.097820 

19.100468 

19.103654 

19.107488 

19.112102 

19.117653 

19.124328 

19.132355 

19.142001 

19.153591 

19.167505 

19.184202 

19.204223 

19.228208 


MAXIMUM 

LOOP 

MAGNITUDE 

(dB) 

39.084840 

39.084838 

39.084835 

39.084832 

39.084829 

39.084824 

39.084819 

39.084812 

39.084805 

39.084795 

39.084784 

39.084771 

39.084755 

39.084735 

39.084711 

39.084683 

39.084649 

39.084608 

39.084558 

39.084498 

39.084426 

39.084339 

39.084235 

39.084108 

39.083957 


MINIMUM 
LOOP 
PHASE 
(degrees ) 

-0.163975 

-0.179963 

-0.197509 

-0.216759 

-0 . 237896 

-0.261099 

-0.286548 

-0.314490 

-0.345153 

-0.378802 

-0.415746 

-0.456281 

-0.500769 

-0.549587 

-0.603176 

-0.661981 

-0.726525 

-0 . 797368 

-0.875115 

-0.960442 

-1.054101 

-1.156881 

-1.269685 

-1.393499 

-1.529391 


MAXIMUM 
LOOP 
PHASE 
( degrees ) 

0.824225 

0.904553 

0 . 992696 

1.089393 

1.195542 

1.312042 

1.439788 

1 . 580001 

1.733815 

1.902522 

2.087651 

2 . 290643 

2.513245 

2.757274 

3.024847 

3.318047 

3.639J25 

3.991236 

4.376505 

4.798082 

5.259197 

5.763069 

6.313262 

6.913456 

7 . 567354 


Table A.l (continued on the next page) 
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FREQUENCY 

( rad /sec ) 

MINIMUM 

LOOP 

MAGNITUDE 

(dB) 

MAXIMUM 

LOOP 

MAGNITUDE 

(dB) 

MINIMUM 

LOOP 

PHASE 

(degrees) 

MAXIMUM 
LOOP 
PHASE 
( degrees ) 

0.102353 

19.256911 

39.083773 

-1.678545 

8.278736 

0.112332 

19.291214 

39.083553 

-1.842244 

9.051234 

0.123285 

19.332152 

39.083286 

-2.021935 

9.888473 

0.135305 

19.380914 

39 . 082966 

-2.219151 

10.793459 

0.148497 

19.438877 

39.082578 

-2.435623 

11.768852 

0.162975 

19 . 507603 

39.082112 

-2.673231 

12.816407 

0.178865 

19 . 588857 

39.081548 

-2.934059 

13.936818 

0.196304 

19.684595 

39.080868 

-3.220372 

15.129205 

0.215443 

19.796956 

39.080046 

-3.534675 

16.390829 

0.236449 

19.928244 

39.079053 

-3.879741 

17.716690 

0.259502 

20.080848 

39.077853 

-4.258567 

19.098797 

0.284804 

20.257226 

39.076400 

-4.674528 

20.526196 

0.312572 

20.459752 

39.074640 

-5.131258 

21.984107 

0.343047 

20.690656 

39.072506 

-5.632812 

23.454113 

0.376494 

20.951883 

39.069917 

-6.183664 

24.913982 

0.413201 

21.244936 

39.066769 

-6.788704 

26 . 337709 

0.453488 

21 . 570794 

39.062936 

-7.453392 

27.696220 

0.497702 

21.929745 

39.058261 

-8.183682 

28.957730 

0 . 546228 

22.321358 

39.052548 

-8.986220 

30.088871 

0 . 599484 

22.744382 

39.045550 

-9.868264 

31.055427 

0.657933 

23.196802 

39.036958 

-10.837894 

31.823537 

0.722081 

23.675826 

39.026379 

-11.903977 

32.360490 

0.792483 

24.177970 

39.013321 

-13.076258 

32.635570 

0.869749 

24.699143 

38.997160 

-14.365425 

32.620679 

0.954548 

25.234739 

38.977105 

-15.783119 

32.2907 W 

1.047616 

25.779736 

38.952163 

-17.341960 

31.623967 

1.149757 

26.328751 

38.921086 

-19.055381 

30.601996 

1.261857 

26.876152 

38.882319 

-20 . 937602 

29.209930 

1.384886 

27.416084 

38.833943 

-23.003270 

27.436410 

1.519911 

27.942551 

38.773620 

-25.267141 

25.273592 

1.668101 

28.449450 

38.698546 

-27 . 743470 

22.717262 

1.830738 

28.930629 

38.605419 

-30 . 445322 

19.766892 

2.009233 

29.379970 

38 . 490434 

-33.383828 

16.425585 

2.205131 

29.791435 

38.349321 

-36.567155 

12.700191 


Table A.l 

(continued on 

the next page) 
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FREQUENCY 

( rad /sec ) 

MINIMUM 

LOOP 

MAGNITUDE 

(dB) 

MAXIMUM 

LOOP 

MAGNITUDE 

(dB) 

MINIMUM 
LOOP 
PHASE 
( degrees ) 

MAXIMUM 
LOOP 
PHASE 
(degrees ) 

2.420128 

30.159146 

38.177417 

-39.999667 

8.601140 

2.65^088 

30.477464 

37.969784 

-43.681168 

4.142189 

2.915053 

30.741031 

37.721377 

-47.606250 

-0.659768 

3.199267 

30 . 944820 

37.427222 

-51.764142 

-5.785181 ' 

3.511192 

31.084155 

37.082611 

-56.138795 

-11.212148 

3.853529 

31.154703 

36.683276 

-60.709375 

-16.916702 

4.229243 

31.152469 

36.225534 

-65.451089 

-22.873000 

4.641589 

31.073781 

35.706374 

-70.336157 

-29.053311 

5.094138 

30.915293 

35.123510 

-75.334820 

-35.427825 

5.590810 

30.674022 

34.475372 

-80.416297 

-41.964428 

6.135907 

30.347414 

33.761079 

-85.549564 

-48.628458 

6.734151 

29.933448 

32.980390 

-90.703934 

-55.382622 

7.390722 

29.430767 

32.133662 

-95.849459 

-62.187130 

8.111308 

28.838810 

31.221804 

-100.957277 

-69 . 000273 

8.902151 

28.157936 

30.246255 

-105.999868 

-75.779242 

9.770100 

27.389503 

29.208968 

-110.951368 

-82.481315 

10.722672 

26.535894 

28.112381 

-115.787971 

-89.065255 

11.768120 

25.600473 

26.959384 

-120.488390 

-95.492766 

12.915497 

24.587485 

25.753273 

-125.034287 

-101.729764 

14.174742 

23.501894 

24.497668 

-129.410736 

-107.747508 

15.556761 

22.349192 

23.196425 

-133.606516 

-113.523305 

17.073526 

21.135182 

21.853524 

-137.614287 

-119.040863 

18.738174 

19.865778 

20. '72970 

-141.430565 

-124.290248 

20.565123 

18.546821 

19.058683 

-145.055550 

-129.267546 

22.570197 

17.183921 

17.614416 

-148.492813 

-133.974281 

24.770764 

15.782340 

16.143675 

-151.748884 

-138.416685 

27.185882 

14.346915 

14.649679 

-154.832787 

-142.604907 

29.836472 

12.882004 

13.135319 

-157.755579 

-146.552232 

32.745492 

11.391467 

11.603144 

-160.529898 

-150.274341 

35.938137 

9.878666 

10.055364 

-163.169557 

-153.788663 

39.442061 

8.346479 

8.493845 

-165.689204 

-157.113821 

43.287613 

6.797318 

6.920130 

-168.104029 

-160.269180 

47.508102 

5.233163 

5.335448 

-170.429525 

-163.274477 


Table A.l 

(continued on 

the next page) 
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FREQUENCY 

( rad /sec ) 

52.140083 
57.223677 
62.802914 
68.926121 
75.646333 
83.021757 
91.116276 
100 . 000000 


MINIMUM 

LOOP 

MAGNITUDE 

(dB) 

3.655588 
' 2.065784 
0.464590 
-1.147495 
-2.770294 
-4.403941 
-6.048876 
-7.705841 


MAXIMUM 

LOOP 

MAGNITUDE 

(dB) 

3.740731 
2 . 136628 
0.523514 
-1.098501 
-2.729566 
-4.370092 
-6.020750 
-7.682473 


MINIMUM 
LOOP 
PHASE 
(degrees ) 

-172.681311 

-174.874974 

-177.025960 

-179.149461 

-181.260324 

-183.372939 

-185.501114 

-187.657931 


MAXIMUM 
LOOP 
PHASE 
( degrees ) 

-166.149540 

-168.914069 

-171.587462 

-174.188680 

-176.736129 

-179.247538 

-181.739836 

-184.229000 


Table A.l : Loop gain extrema at all the trial frequencies, 
for the G(s) constituting a solution for the 
design problem in example 5.3, chapter 5. 


Observe that the fuzzy open loop gain can be perceived to 
cross the negative real axis at around 75 rad/sec with a maximum 
loop gain of less than 1.0. This would mean that the open loop 
gain polar plot does not encircle the critical -l+j0 point. Since 
the given system is minimum phase, the zero encirclement count 
would imply a stable closed loop system. This fact is confirmed 
by the set of stable (i.e. bounded) step responses obtained (see 
figure t. 11, chapter 5) for a randomly selected set of values of 
the uncertain parameters "a" and "b" . 



